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Abstract. 

The MINOS experiment, operating in the NuMI beam since 2005, has provided the most 
precise measurement of the atmospheric mass splitting |Am§ 2 |, and the recent combination 
of the i/ M , and atmospheric neutrino samples has provided some evidence of non-maximal 
mixing, and hints about the neutrino mass hierarchy and the @23 octant. Construction of the 
NOi/A experiment, situated off-axis in the upgraded NuMI beam, is almost complete. Over the 
coming years it will have significant power to probe the questions of the mass hierarchy, #23 
octant, and the possibility of CV violation in the lepton sector. This paper gives an overview 
of the results from MINOS, and of the sensitivity of the NOi'A experiment. 


1. Introduction 

The MINOS and NOi'A experiments are both situated on the NuMI neutrino beam-line. 120GeV 
protons from the Fermilab Main Injector strike a graphite target, the resulting pions are focussed 
by two magnetic horns, and decay to produce a beam of primarily muon neutrinos. Observing 
the energy-dependent disappearance of these muon neutrinos, and subdominant appearance 
of electron neutrinos provides information about the neutrino oscillation parameters. Both 
experiments make use of a two-detector design, in which observations by a Near Detector allow 
measurements of the beam flux, cross-sections and detector effects, while a larger, but otherwise 
similarly constructed, Far Detector observes the change in the beam composition over a much 
longer baseline. 

2. MINOS 

The MINOS detectors are magnetized steel-scintillator tracking sampling calorimeters. They 
consist of alternating planes of 2.54cm thick steel and 4.lcmx lcm solid scintillator strips. Each 
scintillator strip contains an embedded wavelength-shifting fibre, read out at both ends by 
Hamamatsu multi-pixel PMTs. Adjacent scintillator planes are oriented at 90° to each other, 
allowing for 3D event reconstruction. The detectors are magnetized by a coil to a field of around 
IT, allowing for charge identification of muons from curvature. 

The MINOS Far Detector is located underground in the Soudan mine, Minnesota, 735km from 
the NuMI target. It has a total mass of 5.4kton, and has been operating since 2003. The Near 
Detector is located underground at Fermilab, about 1km downstream from the target. It has a 
mass of about lkton, and was completed in 2005. Since 2005, the NuMI beam-line has delivered 
10.7 X 10 20 protons-on-target (POT) to MINOS in the “low energy” neutrino configuration, and 
3.36 x 10 20 POT in antineutrino configuration, plus smaller samples in special configurations. In 
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Figure 1. Energy spectra 
from the MINOS neutrino 
disappearance analyses. The 
black points are the observed 
event counts, with statisti¬ 
cal errors. The grey curves 
show the expectation in the 
absence of oscillations, and 
the red curves show the best 
oscillation fit. Agreement 
with the oscillation hypoth¬ 
esis is very good for both 
neutrino mode (left) and an¬ 
tineutrino mode (right). The 
lower panels display ratios to 
the unoscillated prediction. 


addition, MINOS has collected an exposure of 37.8 kton years of atmospheric neutrinos. In 2013, 
following a shutdown, the NuMI beam returned, delivering neutrinos to NOz^A (and MINOS+ 
and MINERVA) in the new “medium energy” configuration, optimized to give the greatest flux 
around oscillation maximum (~ 2GeV) given NONA’S off-axis location. 

Extrapolating from the observed neutrino spectrum in the Near Detector, in the absence of 
oscillations, MINOS would expect to select 3201 interactions in the Far Detector, and 363 
interactions in the antineutrino running. The observed counts are 2579 and 312 respectively. 
Figure [T] shows the spectral shapes of these deficits, compared to the best fit from the neutrino 
oscillation hypothesis. The three-flavour oscillation paradigm fits the data well: 18% of Monte 
Carlo pseudo-experiments have a worse fit [Tj. 

The inclusion of atmospheric neutrino samples into the analysis provides additional 
sensitivity, particularly to the mass hierarchy. The unoscillated expectation is for 1100 events, 
905 are observed. Figure [2] shows the spectra of the contained muon-selected events as a function 
of zenith angle. Partially-contained and showering events are also included in the analysis as 
additional samples. 

Selecting electron neutrino candidates, the expectation for #13 = 0 is 69.1 events of 
background in neutrino mode, and 10.5 in antineutrinos. The counts observed in data are 
88 and 12. Figure [3] shows these excesses in bins of energy and PID value. The hypothesis that 
$13 = 0 is rejected at the 96% confidence level [2j. 

Combining the information provided by the atmospheric, and v e samples, plus additional 
constraints on $13 from reactor experiments provides some sensitivity to the neutrino hierarchy 
and 023 octant. The solar parameters (AtttJi, $ 12 ) are held fixed in this analysis. $13 is treated 
as a nuisance parameter, constrained by external reactor results. Scp, 023, and Am | 2 are 
unconstrained and determined by the joint best fit to all the spectra. Major sources of systematic 
error are included as nuisance parameters. 

Figure [4] shows the confidence contours obtained in (Am| 2 , sin 2 623 ) space. The best fit 
value of sin 2 023 is 0.41, and maximal mixing is disfavoured at 79% confidence. The data also 
weakly favour the inverted hierarchy and the lower O 23 octant. Figure [5] shows a different 
projection of this information, with 5cp information included. The most disfavoured scenario is 
the combination of normal hierarchy and upper O 23 octant which is rejected at the 81% confidence 































Figure 2. Spectra for the MINOS atmospheric neutrino sample. The points show the observed 
event counts as a function of reconstructed zenith angle. The grey curves show the unoscillated 
predictions, and the red and blue dashed curves are the best oscillation fits assuming normal 
or inverted hierarchy respectively. The filled blue histograms indicate the background due to 
cosmic-ray muons. 


Table 1. Summary of MINOS combined oscillation fit results for normal and inverted hierarchy. 


Parameter 

Best fit 

Confidence limits (NH) 

Best fit 

Confidence limits (IH) 

|Ara 2 2 |/10 _3 eV 2 

2.37 

2.28 - 2.46 (68% C.L.) 

2.41 

2.32 - 2.53 (68% C.L.) 

sin 2 $23 

0.41 

0.35 - 0.65 (90% C.L.) 

0.41 

0.34 - 0.67 (90% C.L.) 


Preference for inverted hierarchy: — 2AlogL=0.23 
Preference for lower octant: — 2AlogL=0.09 
Exclusion of maximal mixing: —2AlogL=1.54 (—>• 79% C.L.) 


level for all values of 5 cp- This measurement is also the most precise determination of |Am| 2 |. 
The results are summarized numerically in Table |T] [3]. 

If some component of the disappearing MINOS flux were oscillating to sterile neutrinos, one 
would expect to also see a deficit of neutral current events, in addition to the deficit seen in 
charged currents. In fact, a search for this effect found a small excess of neutral current 
interactions at the Far Detector as compared to predictions. Interpreting in terms of oscillations 
around A m? ~ 0.5eV 2 leads to a limit sin 2 2 0^ e < 7.1 x 10~ 3 at the 90% confidence level. 

The MINOS detectors continue to operate in the NOi'A-era beam. The higher energy peak 
of the neutrino flux reduces the statistics acquired in the oscillation region , but the statistics 
at higher energies are greatly increased. MINOS+ expects to collect around 4,000 charged 
current events per year. This high-statistics sample will allow tests of the standard three-flavour 
oscillation paradigm, including searches for sterile neutrinos, and sensitivity to signatures of 
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Figure 3. Spectra of ^-selected 
events in MINOS. The points are 
the number of events selected as 
a function of reconstructed neu¬ 
trino energy. The red curve 
shows the prediction for $13 = 0, 
and the purple region shows the 
excess for the best-fit value of 
013. The samples are divided by 
PID value (top to bottom) and 
between neutrino (left) and an¬ 
tineutrino (right) running. The 
data is divided into bins of PID 
value, from the lowest selected 
PIDs (top) to the most i' e -like 
events (bottom). 88 neutrino¬ 
mode events and 12 antineutrino¬ 
mode events are selected, com¬ 
pared to an expectation for $13 = 
0 of 69.1 and 10.5. The expecta¬ 
tions for 0i3 = 0.1 are 95.1 and 
13.6. 
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exotic physics such as large extra dimensions. 

3. NOz^A 

The NOzxA. detectors are fine-grained, low-Z, highly active tracking calorimeters. Cells of 
extruded PVC, approximately 6cm x 4cm are formed into planes, with the orientation of 
successive planes alternating as in MINOS. They are filled with a mineral oil/liquid scintillator 
mix, which comprises 64% of the detector by mass. Each cell contains a looped wavelength- 
shifting fibre, read out at one end by one of the 32 pixels of a Hamamatsu avalanche photodiode 
(APD), cooled to -15°C @3- 

The Far Detector is located on the surface at Ash River, Minnesota, 810km from the NuMI 
target. The detector consists of 344,000 channels, approximately 16 x 16 x 60m, for a total mass 
of 14kton. Assembly and placement of the plastic blocks, scintillator filling, and electronics 
installation proceed in parallel. All blocks are now in place, and oil filling is almost complete. 
Over 2/3 of the detector is instrumented with APDs and being read out. The Near Detector 
is located underground at Fermilab, approximately lkm from the target. It comprises 18,000 
channels, with a mass of 0.3kton. All blocks are in place, and oil filling and electronics outfitting 
are in progress. The NuMI beam is being upgraded from a maximum capacity of 350kW up to 
700kW. This will be achieved by converting the Recycler from running antiprotons to protons, 
and shortening the cycle in the Main Injector from 2.2s to 1.33s. The beam has been running 
stably since returning from shutdown in September 2013, and power is gradually being increased. 
The maximum power achievable will be limited to 500kW until upgrading of the Booster RF 




























































Figure 4. Combined results from the MINOS v e , and atmospheric analyses. The star 
shows the best fit values of sin 2 023 and Am| 2 - This point weakly favours non-maximal mixing, 
the lower octant, and inverted hierarchy. The red and blue contours indicate 68% and 90% 
confidence levels. The panels on the right display the profile likelihood projection onto each of 
the two parameters. 


cavities is complete. Thus far, in excess of 1.5 x 10 20 protons have been delivered. 

It is important to validate the detectors’ ability to pick neutrino interactions out of the cosmic- 
ray background, and to verify the synchronization of the Far Detector to the NuMI spill times. 
Analysis of the Near Detector prototype on the surface at Fermilab (NDOS) shows a clear peak 
in the neutrino candidate timing spectrum over the cosmic-ray background. Propagating this 
timing information to the Far Detector gives a lO^us window in which neutrinos are expected to 
be found. A combination of automated analysis and hand-scanning of events in a 500/us window 
around this point (plus out-of-time spills for background estimation) found four candidate events 
within the timing window, on a background of 0.05. The times of the events relative to the beam 
window were kept secret until after the choice of these events had been finalized. Figure [7] shows 
the NDOS timing peak compared to the times of the selected neutrino candidates, and Figure 
[8] shows one of the Far Detector events in the timing window, a multi-pronged interaction, with 
reconstructed neutrino direction coming from Fermilab. 

The sensitivity projections presented below assume 18 x 10 20 POT of data taken in neutrino 
mode plus 18 x 10 20 POT taken with antineutrinos. This corresponds to a nominal 6 x 10 20 POT 
per year and a run plan with three years in each mode. The choice of what exposure to take in 
each beam configuration is flexible, and may be adjusted in response to changing circumstances. 
The oscillation parameters assumed are sin 2 2013 = 0.095 sin 2 2023 = 0.95 or 1.0 depending on 
context. We take AmJ 2 = 2.35 x 10 _3 eV 2 , and the other parameters from the latest Particle 
Data Group averages [5]. 

Current global constraints on the atmospheric mixing angle are around sin 2 2023 > 0.9, 














Figure 5. Log-likelihood differ¬ 
ences from the best fit oscillation 
parameters for different values of 
8 cPi hierarchy and 023 octant for 
the combination of MINOS analy¬ 
ses. The region normal hierarchy, 
upper octant, 5cp ~ 4r is dis¬ 
favoured at greater than 90% con¬ 
fidence. 


Figure 6. Sensitivity of NOi'A for a nominal run of 
three years neutrino plus three antineutrino (18 x 10 20 
POT each). The true oscillation values are chosen to be 
close to the MINOS best-fit. The shape of the curves 
is very similar to the MINOS results in Figure 0 but 
the expected sensitivity is much higher. For example, 
the most disfavoured scenario, normal hierarchy, upper 
octant, would be excluded at 5<r for all values of 5c p- 
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Figure 7. The timing peak formed by neutrinos 
from the NuMI beam in the NDOS detector 
(filled histogram) with the times of selected 
interactions in the NOi'A Far Detector marked 
(circles). The Far Detector events were found by 
a combination of hand-scanning and automated 
analysis. The estimated cosmic-ray background 
in the 10^s Far Detector spill window is 0.05 
events. 


with hints of non-maximal mixing from MINOS. The NOz^A experiment will have significantly 
enhanced power to distinguish between maximal and non-maximal mixing. The u tl charged 
current analysis being developed separates the data into three samples: quasi-elastic candidates, 
which are very pure with good energy reconstruction; the remainder of contained candidates; 
and uncontained events, where the muon exits through the back or side of the detector, a sample 
with substantial statistics, but poor energy resolution. Predicted event counts are given in Table 
[21 The contours that would be obtained with increasing exposure are shown in Figure [9j for two 
possible values of $23- With three years of neutrino running, and three years of antineutrinos, 
this analysis should obtain percent-level uncertainty on the atmospheric mixing parameters. If 
sin 2 023 = 0.95, maximal mixing could be excluded at 90% confidence with one year of running. 

Table [2] also shows the representative event counts for the sample of electron neutrino 
candidates. As with MINOS, sensitivity to the mass hierarchy, 023 octant, and 8 cp comes from 
combining the disappearance and appearance analyses, plus constraints on 0i3 from reactor 
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Figure 8. Event display of the first neu¬ 
trino candidate event selected in the NOi'A 
Far Detector. Each square represents en¬ 
ergy deposited and detected in a single 
cell. The neutrino beam enters from the 
left. The top panel shows the informa¬ 
tion from vertically-oriented cells, with hor¬ 
izontal cells shown in the bottom panel. 
The overlaid coloured lines are the recon¬ 
structed paths of particles, and the cross 
marks the reconstructed position of the in¬ 
teraction vertex. 
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Table 2. Representative event counts for NOz^A analyses. The columns represent 18 x 10 20 
POT neutrino-mode and antineutrino-mode running. The analysis divides the dataset into 
three samples, candidate quasielastic events, non-quasielastic events, and events where the muon 
is not fully contained within the detector. The counts are constrained to the range 0-5GeV, the 
region of sensitivity to oscillations driven by the atmospheric mass splitting. 
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NC 

19 
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vn CC 
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Beam v e 

8 

5 

Tot bkg 

32 

15 
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68 

32 
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CC 

NC 
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NC 

Quasielastic 

82 

< 1 

49 

< 1 

non-QE 

168 

14 

78 

6 

Uncontained 

233 

6 

134 
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experiments. Figure [TO] shows the sensitivity of this analysis to these three parameters, as 
a function of the true values of 5c p and the mass hierarchy. For the mass hierarchy, there 
are degenerate regions where it is not possible to disentangle the matter effects, which give 
information about the mass hierarchy, from true CP-violation driven by 5cp- But for favourable 
combinations 2 — 3cr evidence for the mass hierarchy can be obtained. Determination of the #23 
octant is much less sensitive to 5cp and can be achieved at better than 2cr for all values of 5cp 
if sin 2 2023 = 0.95. The lower octant is assumed in order to match the best fit from MINOS, 
performance for the upper octant is slightly better. Figure [6] shows predicted Ay; 2 slices as 
a function of 5cp, to be compared with the analogous Figure [5] from MINOS data. The true 
parameters are chosen to match the MINOS best fit (inverted hierarchy, lower octant, 5cp ~ §)• 













Figure 9. Projected 90% confidence 
contours in the (Am| 2 , sin 2 # 23 ) space for 
two, six, or ten years of running (assuming 
6 x 10 20 POT per year), for two values 
of 6 * 23 . In the case of sin 2 2023 = 0.95, 
maximal mixing can be excluded at 90% 
confidence in under two years of running. 
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Figure 10. Sensitivity of NOi'A to resolve the 
mass hierarchy (top left), discover CP violation 
(above), and resolve the 023 octant (left) as 
a function of the true value of Scp- The 
solid blue curves are in the case that the 
true mass hierarchy is normal, the dashed red 
curves are for the inverted case. In each case, 
sin 2 2013 = 0.095 is assumed. For hierarchy and 
CP violation sin 2 2023 = 1 is assumed, the octant 
sensitivity assumes sin 2 2023 = 0.95, first octant. 


If the true parameters are indeed in this region, NOz^A would reject parts of phase space (in 
particular the combination of normal hierarchy and upper octant) at very high significance. 

4. Conclusion 

The MINOS experiment has collected a large sample of neutrino interactions in over six years of 
running, in both neutrino and antineutrino mode. A combined fit to this data gives the world’s 
most precise measurement of |Am| 2 |, and hints for the neutrino mass hierarchy and 623 octant. 
These detectors will continue taking data in the NOi^A-era beam as MINOS+. 

Construction of the NOz/A experiment is progressing well. First neutrinos have been found 
in the Far Detector, and projections show good sensitivity to the mass hierarchy and 623 octant, 
in addition to improved measurements of the atmospheric mixing parameters. Results from the 
next few years, in combination with other experiments, should provide significant information 






































































on these questions for the first time. 
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